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The diatomic linear chain of masses coupled by harmonic 
springs is a textboook model for vibrational normal modes 
(phonons) in crystals. In addition to propagating acoustic 
and optic branches, this model is known to support a "gap 
mode" localized at the surface, provided the atom at the sur- 
face has light rather than heavy mass. An elementary argu- 
ment is given which explains this mode and provides values 
for the frequency and localization length. By reinterpreting 
this mode in different ways, we obtain the frequency and lo- 
calization lengths for three other interesting modes: (1) the 
surface vibrational mode of a light mass impurity at the sur- 
face of a monatomic chain; (2) the localized vibrational mode 
of a stacking fault in a diatomic chain; and (3) the localized 
vibrational mode of a light mass impurity in a monatomic 
chain. 



I. INTRODUCTION 

Unlike molecules which have discrete vibrational fre- 
quencies, crystals have a continuous spectrum of vibra- 
tions which can propagate as travelling waves |l] |3| . This 
fact causes crystals to be much better heat conductors 
than glasses or liquids. Sometimes the spectrum is in- 
terrupted by gaps where no propagating normal modes 
occur. Other interesting behavior happens at frequen- 
cies inside the gap, such as localized (non-propagating) 
normal modes associated with defects and surfaces. The 
text by Ziman Q has a good discussion. A visualization 
of surface modes on the (100) surface of Cu is on the 
website of Ch. W611, Ruhr-Universitat Bochum j5j. The 
present paper shows how this happens for some partic- 
ular cases of one-dimensional crystals, or linear chains 
of atoms. Our treatment uses only classical mechanics, 
and gives properties (frequency and displacement pat- 
tern) rigorously by pictorial arguments with no higher 
algebra. 

The surface phonon provides the simplest example 
of wave localization, an effect which occurs in many 
branches of physics. Analogous phenomena are found in 
the quantum treatment of electrons in single-particle ap- 



proximation [M 2 1 , and in the new field of "photonic band- 
gap systems" [B]. This paper reports a simple way of 
understanding the surface phonon on the diatomic linear 
chain. The model is then extended and reinterpreted to 
give simple explanations of some other localized modes. 



II. DIATOMIC MOLECULE 

A diatomic molecule has a single vibrational "normal 
mode." Even though the restoring force of atom 1 on 
atom 2 has in reality a complicated quantum-mechanical 
origin, for small displacements away from equilibrium, 
it can always be well approximated by a spring obeying 
Hooke's law with a spring constant K. Using standard 
physics of the two-body problem jjj, if the two atoms 
have masses Mjj and (H and L are for heavy and 
light), the squared oscillation frequency to 2 is K/M Ye< i 
where M red is the "reduced mass" M L M H /(M L + M H ). 



III. PERFECT INFINITE CHAINS 

A crystalline solid is a very large molecule, with a con- 
tinuous spectrum (or band) of vibrational frequencies. 
Solids can also be modelled by masses connected to each 
other by springs. A one-dimensional chain of masses is 
often studied, not because it is found in nature, but be- 
cause the mathematics is simple and can be generalized 
to more realistic three-dimensional arrangements. For a 
large enough collection of atoms, most of the vibrational 
normal modes are classified as "bulk" normal modes, 
which means they are essentially identical to those of 
a hypothetical infinite sample with no boundaries. Each 
"bulk" normal mode has a pattern of atomic displace- 
ments which extends throughout the system. Similar to 
the normal modes of a vibrating string, these are sine 
and cosine standing waves. Alternately, one can use lin- 
ear combinations of sines and cosines to give an equiva- 
lent basis of left and right-going travelling waves. For the 
case of all masses equal to Mo, the Pth atom (located at 
Re = £a) has a displacement Aam(kRe —tOkt) in a right- 
going travelling wave. The corresponding squared fre- 
quency is (4K/M )sm 2 (ka/2). There are as many such 
solutions (N) as there are atoms in the chain, namely 
solutions for each k in the range (— 7r/a, ir/a). This is 
derived in many texts fl]-|||§ ■ For N — > oo the spectrum 
is continuous between the minimum squared frequency 
of zero and the maximum of w^ax = ^K/Mq. A partic- 
ularly original discussion is given by Martinez 0. 

The vibrational spectrum of a real material sometimes 
has a gap, an interval of frequencies where there are no 
travelling wave solutions. A simple model illustrating 
this is the "diatomic chain," an infinite chain of alter- 
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nating masses Ml, Mh- The algebra, which is more 
complicated than the monatomic chain, is also given in 
texts [0,| . The dispersion curve for lo\ is given in Fig. 
0. There are now two "branches," labelled acoustic and 
optic, and a gap. Exactly in the middle of the gap, the 
surface may induce a localized vibrational normal mode, 
with amplitude which falls exponentially (cx exp(— R/£)) 
with distance R into the bulk. 

Before discussing this, we sharpen our understanding 
with a quantitative interpretation of the four special bulk 
modes indicated by circles in Fig. 0. The frequencies 
of these special modes can be understood without the 
algebra needed to find the frequencies of the modes at 
general A:- vectors. 
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frequency (in 



loq — K/M) versus wavevector k for the diatomic chain 



units 
The 

gap is proportional to (Mh — Ml ) /(Mh + Ml ) , and is drawn 
for the case Mh = 2Ml- The distance between atoms is a. 
The dashed line shows the position of the surface mode. 



mode b: This has oppositely directed displacements for 
adjacent atoms. Each unit cell of the crystal has 
the same displacment pattern. Therefore the wave- 
length is infinite and the wavevector is zero. The 
displacements in mode b arc such that u L (the 
displacement of the light atom) is proportional to 
Mh, and similarly uh is proportional to Ml- Thus 
the center of mass of each unit cell is fixed. The 
mode is almost the same as in a diatomic molecule, 
except each atom has two springs attached, one 
stretched and the other compressed by the same 
amount. Therefore, when released from rest, each 
pair of atoms oscillates with fixed center of mass 
but with twice the restoring force of an isolated di- 
atomic molecule, i.e. cu 2 = 2K/M ro d- This is the 
highest frequency normal mode in the spectrum. 

mode c: This has light atoms stationary and heavy 
atoms moving in an alternating pattern. The light 
atoms feel equal and opposite forces which cancel, 
while the heavy atoms feel repulsive and attractive 
forces which add. This pattern also oscillates in 
time, with squared frequency lu 2 — 2K /Mh ■ 

mode d: This is the same as mode c except heavy and 
light atoms are interchanged, making the squared 
frequency equal to 2K/Ml- Modes c and d have 
wavelength 4a and wavevector it /2a. All other nor- 
mal modes of the infinite crystal are more com- 
plicated and have frequencies which lie on smooth 
curves connecting these four modes. 
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IV. SPECIAL BULK MODES 



The four special modes circled in Fig. 0have the sim- 
ple vibrational patterns shown in Fig. First, why 
are these patterns "normal modes" ? If we take as initial 
conditions, the velocities of all atoms to be zero and the 
positions to be as shown in the figures, then Newton's 
laws have simple, and perhaps even obvious solutions: 
the pattern is preserved, and oscillates in time as cos(wi) 
for some special choice of ui. This is the definition of 
a normal mode. Second, what is the corresponding fre- 
quency of oscillation? This can be answered by careful 
consideration of forces and masses. 

mode a: This is the simplest mode with all atoms hav- 
ing the same displacement. This has infinite wave- 
length (zero wavevector), no stretch of any spring, 
and therefore zero restoring force and w = 0. 
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FIG. 2. The diatomic chain, the special bulk modes a, b 
with wavelength A = oo and wavevector k — 0, the special 
bulk modes c, d with wavelength A = 4a and k = it /2a (all 
shown circled in Fig. Q), and the gap mode e confined at the 
surface when a light atom terminates the chain. In pictures 
a d, the chain is infinite; by contrast, in picture e, the light 
atom at the left is the surface atom, with no spring acting 
to its left. The brackets indicate pairs of atoms which move 
without altering the length of the bond between them. 
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V. SURFACE MODE IN THE GAP 

Modes which are confined to the surface region nor- 
mally must have frequencies which lie outside the "bulk" 
bands. Discussions of such modes are given in texts on 
surface physics ]ic|-|l^] and measurements are cataloged 
by Kress and de Wette jlB). We have discovered a very 
simple explanation of the fact |l6| that a "gap mode" 
confined to the surface occurs in the diatomic chain if 
the endmost atom is a light atom. 

Consider mode e, which like mode b has pairs of atoms 
vibrating with fixed center of mass. However, adjacent 
pairs vibrate in such a way that the connecting spring 
is not stretched. Thus each pair experiences no force 
from any other atom and is decoupled from the rest of 
the chain. The resulting decoupled pairs oscillate with 
u! 2 = K/M Ic d as for isolated diatomic molecules. Since 
all pairs have the same frequency this is a stable nor- 
mal mode. The frequency lies exactly in the middle of 
the gap of the squared frequency spectrum (K/M Te A = 
(l/2)(2K/M L + 2K/M H )). In order to be decoupled, 
the heavy atom of a given pair, and the adjacent light 
atom of the next pair deeper into the bulk, must have 
the same displacement, smaller by Ml/Mh (and with 
opposite sign) than the displacement of the previous light 
atom closer to the surface, in order to conserve center of 
mass position. Since adjacent pairs have displacement 
ratios —Ml/Mh, the n-th pair has amplitude propor- 
tional to (—Ml/Mh) 11 = (-l)"exp(-nln(M ff /M L )). 
This is an exponential decay exp(— 2na/£) with decay 
length £ = 2a/ \n(Mn /Ml). If the surface atom had 
been a heavy atom, this mode would have been exponen- 
tially growing rather than decaying, which is not allowed 
for a normal mode. Mode e was first found by Wallis @ 
in an elegant calculation of the spectrum of finite chains. 
Our simple argument is not (to our knowledge) in the 
literature. A "standard" derivation is given in the text 
by Cottam and Tilley Q. 

VI. SURFACE MODES OF 
THREE-DIMENSIONAL CRYSTALS 

Mode e is directly related to a branch of surface nor- 
mal modes of higher-dimensional diatomic crystals. A 
two-dimensional version is shown in Fig. |^. Various types 
of surfaces are possible for such crystals. If cut perpendic- 
ular to a conventional x or y axis shown in the figure by 
dashed lines, the surface contains equal numbers of A and 
B ions, and is referred to as "non-polar." By contrast, 
the surface shown is a "polar surface" with a layer of 
light atoms exposed and layers of heavy and light atoms 
alternating underneath. There is a vibrational normal 
mode in which each layer oscillates perpendicular to the 
surface (as indicated by arrows) and which is localized 
at the surface. Of course, in real crystals the forces ex- 
tend beyond first neighbors, so the displacement ratio 



(—Ml/Mh) may not be exactly obeyed and the squared 
frequency may not lie exactly at mid-gap, but the actual 
behavior will mimic reasonably well the idealized one- 
dimensional example of the previous section. 
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FIG. 3. "Polar" surface of a diatomic crystal in two di- 
mensions. This is a two-dimensional analog of the "(111)" 
surface of crystals with the "rocksalt" (NaCl) structure. A 
square conventional unit cell is shown as dashed lines parallel 
to conventional x and y axes. There are two possible ter- 
minations of such crystals, ones which expose planes of light 
atoms as shown in the figure by open circles, and ones which 
expose planes of heavy atoms, shown as filled circles. Since 
the heavy and light atoms have opposite ionic charges, such 
surfaces have a surface dipole and are called polar. 

There is actually not just one mode of this type, but a 
branch of such modes, with displacement patterns sinu- 
soidally modulated along the surface. The one depicted 
in Fig. H has the surface atoms "a" , "b" , "c" , all mov- 
ing in phase, corresponding to an infinite wavelength, 
or zero wavevector, parallel to the surface. The other 
extreme case of modulation is when atoms along the sur- 
face are completely out of phase; when atom "a" moves 
down, atom "b" moves up, and so forth, corresponding 
to a wavelength A = 202a in the plane of the surface. 
Thus we anticipate a branch of surface excitations with 
wavevectors lying in the plane of the surface. In order for 
such a mode to be exponentially localized in the surface 
region, the frequency of oscillation must lie in a gap where 
there are no corresponding bulk normal modes with the 
same components of wavevector in the plane of the sur- 
face. A gap is almost certain to occur for the case of zero 
wavevector, but at increasing wavevectors the gap may 
disappear, and the mode ceases to be localized near the 
surface. 

Dimension two or three also opens new possibilities 
less directly related to one-dimensional models, such as 
surface normal modes with displacements in the plane 
of the surface. Many branches of surface normal modes 
have been seen experimentally by scattering experiments. 
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Unfortunately we have not been able to locate in the 
literature any observation of the simple mode illustrated 
in Fig. ||. This is perhaps because polar surfaces are 
relatively unstable and hard to create and work with. 



VII. IMPURITY ATOM ON THE SURFACE 

Another known result is that a surface mode appears 
above the bulk frequency spectrum for a monatomic 
chain, provided the atom on the surface is lighter than 
the rest by at least a factor of two. This can be proven 
by a reinterpretation of the previous construction. For 
mode e in Fig. g, let the two atoms connected by the 
unstretched spring be reinterpreted as a single atom of 
mass Mo = Mb + Ml- Then the model has new in- 
terior atoms all with mass Mo, but a surface impurity 
atom with mass Mi mp = Ml < 0.5Mo- The surface 
mode e still solves Newton's laws with lo 2 s — K/M Te ^ 
and M ro d = MlMb / '{Ml + My). In terms of the new 
variables M imp and M the reduced mass M rc( j is M = 
M; mp (Mo — Mi mp )/M . The frequency w| lies above the 
top of the bulk band (wmax = 4K/M ) if M > 2M imp , 
and merges into the bulk band for Mo < 2Mj mp . This 
result seems also to have been first discovered by Wallis 
jl7j . A "standard" proof of this result is in the book by 
Dcsjonqueres and Spanjaard [O. 



VIII. LOCALIZED GAP MODE OF A STACKING 
FAULT 

The gap mode e of Fig. |^ generates a corresponding 
mode of a defective bulk crystal, shown in Fig. [|. This 
mode decays exponentially in both directions away from 
the center of symmetry. This center lies in the middle 
of a "stacking fault" where two light-mass atoms have 
been put adjacent to each other. It is a one-dimensional 
version of a planar defect which occurs in real three- 
dimensional crystals. The quantum-mechanical force be- 
tween two light-mass atoms differs from the force which 
binds the atoms of unlike mass. Therefore, we must ex- 
pect that the separation a 1 of the light-mass atoms will 
differ from the equilibrium separation a of unlike atoms, 
and that the force constant K' between these atoms will 
differ from the constant K occuring elsewhere. Notice 
that for the special displacement pattern of Fig. ^, there 
is no force between the adjacent light atoms, so the val- 
ues of a' and K' are irrelevant; the squared frequency of 
the normal mode is exactly the same as the surface mode 
e of Fig. [^, and is pinned at midgap. 
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FIG. 4. Localized normal mode of vibration of the diatomic 
chain with a stacking fault, obtained from mode e of Fig. ^ 
by reflecting the lattice and assigning u(—l) = u(£) 

The stacking fault is a simple example of a "topolog- 
ical defect," that is, a defect which cannot be trans- 
formed away by any local change. As far as we know, 
the mid-gap normal mode of vibration found here for the 
stacking fault has not previously been discussed in the 
literature. However, a close analog is the "topological 
soliton" found at mid-gap in the electronic spectrum of 
the "Su-Schrieffer-Heeger" model [H for polyacetylene 
with a topological defect in the pattern of dimerization 
of carbon-carbon bonds along the chain. 



IX. LOCALIZED VIBRATION OF A LIGHT 
MASS IMPURITY IN A MONATOMIC CHAIN 

Suppose an impurity of mass Mi mp < Mo is substituted 
into a monatomic chain of mass Ma with no change in 
force constants. Define the fractional mass deficit to be 
e = (Mo — Mj mp )/Mo > 0. It is known that this system 
supports a localized mode whose frequency "splits off" 
above the frequency Umax of the uppermost bulk mode. 
Specifically, the mode has squared frequency w^ AX /(l — 
e 2 ) and is localized around the impurity with localization 
length a/ ln((l + e)/(l — e)). The earliest presentation of 
this mode known to us is by Montroll and Potts Jl9| . 
The topic of localized modes in solids had been given a 
systematic formulation in three earlier papers by Lifshits, 
available only in Russian A textbook derivation is 
given by Mihaly and Martin [BCj, a nd a nice qualitative 
discussion is given by Harrisonj|22j . 

These results follow rigorously by reinterpretation of 
Fig. |J. Simply regard each pair of co-moving atoms as a 
single atom whose mass is the sum of the two shown in 
the figure. Thus M is M H + M L , M imp is 2M L , and the 
new lattice constant a is twice the previous distance a. 
When the impurity mass is heavier than the host mass, 
there is no longer a split-off bound state, but instead a 
"resonance" within the bulk band. 

In three-dimensional crystals the occurrence of a vi- 
brational bound state requires a minimum mass deficit e 
which is model-dependent, whereas our 1-d example has 
a bound state for arbitrarily small mass deficit. This 
is a classical discrete-system analog of the continuum 
quantum-mechanical theorem that an attractive well al- 
ways has a bound state in a 1-d one-electron problem 
(and also in 2-d) but requires a critical well-depth in 3-d 
p3| . For the impurity on the surface, however, we saw 
that even in Id there is a critical mass deficit of 1/2. 
The quantum analog is that if the well is at the edge of 
a 1-d half space (the other half of space is impenetrable 
because of an infinite potential), then there is a critical 
well-depth, equal to the well-depth at which the second 
bound state appears for the symmetric well in the full 
1-d space. 
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X. SUMMARY 

Two simple surface phonons and two simple bound de- 
fect modes in one-dimensional lattices have been quanti- 
tatively explained by pictorial construction and elemen- 
tary physics of the two-body problem. This is certainly 
not a complete catalog of interesting localized modes, but 
we think that these modes can serve as useful pedagogical 
models for phenomena in several branches in physics. 
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